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Abstract. Close binaries consisting of a main sequence 
star and a white dwarf are considered as candidates for 
Type la supernova progenitors. We present selfconsistent 
calculations of the time dependence of the structure of the 
main sequence star, the mass transfer rate, and the orbit 
by means of a binary stellar evolution program. We obtain, 
for the first time, a complete picture of the time evolution 
of the mass transfer rate in such systems. We find a long 
switch-on phase of the mass transfer, about 10® yr, dur- 
ing which nova outbursts should persist in all systems. 
Furthermore, we find that the white dwarfs can reach the 
Chandrasekhar mass only during the decline phase of the 
mass transfer, which may have consequences for the crit- 
ical accretion rate for stationary nuclear burning on the 
white dwarf surface. In contrast to results based on sim- 
ple estimates of the mass transfer rate in systems of the 
considered type, our results allow for the possibility that 
even systems with rather small initial white dwarf masses 
(~ 0.7 Mq) may produce Type la supernovae, which then 
might originate from very rapidly rotating white dwarfs. 

We present results for two different metallicities, 
Z=0.02 and Z=0.001. We find that for systems with the 
lower metallicity, the mass transfer rates are on average 
five times larger than in comparable system at solar metal- 
licity. This leads to a systematic shift of the supernova la 
progenitor population. Firstly, while for Z=0.02 — for our 
choice of white dwarf wind mass loss and mass accumula- 
tion rate — donor star initial masses in supernova progen- 
itor systems are restricted to the range 1.6 M0...2.3 Mq, 
they are in the interval 1.4 Mq...1.8 Mq at low Z. Sec- 
ondly, the initial white dwarf masses need, on average, to 
be larger by 0.2 Mq at low Z in order to obtain a Chan- 
drasekhar mass white dwarf. This metallicity dependences 
have very little effect on the progenitor life times, but may 
be responsible for a drop of the Type la supernova rate 
for low metallicity, and may introduce a Z-dependencc in 
the properties of supernovae which stem from close main 
sequence star + white dwarf systems. 
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We estimate the X-ray luminosities of the computed 
systems, and investigate their donor star and orbital prop- 
erties. We find the donor stars to be underluminous by up 
to one order of magnitude, and more compact than normal 
main sequence stars. In general, our systems correspond 
well to observed close binary supersoft X-ray sources. We 
further derive the chemical and kinematical properties of 
the stellar remnants of our systems after the explosion of 
the white dwarf, which may serve as a test of the viability 
of the considered Type la supernova scenario. 

Key words: stars: abundances - binaries: close — stars: 
evolution - stars: interiors - supernovae: general - X-rays: 
stars 



1. Introduction 

During the last years, the refinement of supernova obser- 
vations, e.g., the routine detection of supernovae at large 
redshifts, has made them a powerful tool to probe cos- 
mology. It allowed to determine the Hubble constant with 
unprecedented accuracy (Riess et al. 1995, Hamuy et al. 
1996; see also Hoflich & Khokhlov 1996). Even more ex- 
citing, recent results (e.g., Riess et al. 1999, Perlmutter 
et al. 1999) are consistent with a low matter density in 
the Universe and, intriguingly, hints for a positive cos- 
niological constant. These findings are based on empirical 
brightness-decline relations which are calibrated locally. 
This leaves potential systematic effects of supernova la 
properties with redshift as major concern. To this end, it 
would be desirable to obtain an estimate of such effects 
from theoretical models of supernova la progenitor sys- 
tems. 

However, despite considerable efforts during the last 
decades, the exact nature of supernova la progenitors is 
still unclear. On observational and theoretical grounds, it 
is generally agreed that Type la supernovae result from 
the thermonuclear disruption of a CO white dwarf (e.g., 
Woosley & Weaver 1986, Wheeler 1996, Nomoto et al. 
1997, Branch 1998). Since isolated white dwarfs cool, a 
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close binary component which transfers mass to the white 
dwarf is a prerequisite to obtain a Type la supernova. 
Various binary evolution scenarios leading to exploding 
CO white dwarfs have been proposed and investigated, 
but hitherto it is unclear which of them is preferred in 
nature (cf. Branch 1988, Livio 1999). 

In this paper, we study the evolution of close binary 
systems consisting of a carbon-oxygen white dwarf and 
a main sequence star, which was repeatedly proposed as 
promising supernova la progenitor scenario (cf. Nomoto 
& Sugimoto 1977, Li & van den Heuvel 1997, Kato & 
Hachisu 1999, Hachisu et al. 1999). In these systems, the 
carbon-oxygen white dwarfs are the remainders of stars 
with an initial mass below « 10 M© which have lost their 
H/He-rich envelope, with CO cores of « O.6...1.2M0. If 
accretion is sufficiently fast the accreted hydrogen may 
burn to helium and, subsequently, to CO on the surface 
of the white dwarf, and its mass grows close to the Chan- 
drasekhar mass. 

The binary evolution leading to close white 
dwarf + main sequence star systems is not yet well 
understood (cf., Livio 1996). However, we know a large 
number (~ 10^) of close white dwarf -|- main sequence 
star systems as Cataclysmic Variables (Ritter & Kolb 
1998), most of which do not evolve into Type la su- 
pernovae since they undergo nova outbursts which may 
prevent a secular increase of the white dwarf mass (e.g., 
Kovetz & Prialnik 1997). The idea that also the slightly 
more massive systems of the same type studied here 
occur in nature is supported by population synthesis 
studies, which predict their birth rate to be comparable, 
within an order of magnitude, to the observed rate of 
la supernovae (e.g., de Kool & Ritter 1993, Rappaport et 
al, 1994). It is further supported through the discovery of 
the so called supersoft X-ray sources (Greiner et al. 1991, 
Kahabka & van den Heuvel 1997), which may represent 
the observational counterparts of the binary systems 
studied here theoretically. 

We investigate the properties of close main sequence 
star-white dwarf systems at two different metallicitics. As 
we derive the detailed time-dependence of the accretion 
rate, our work is relevant for the understanding of indi- 
vidual supernovae and supersoft X-ray binaries, for the 
change of their average properties with metallicity, and for 
the dependence of the rate of la supernovae with metal- 
licity. We introduce our computational method in Sect. 2, 
and present our results for the mass transfer rate and re- 
sulting maximum white dwarf masses in Sect. 3. In Sect. 4, 
we discuss the evolution of the white dwarf spin, of the bi- 
nary orbit and of the main sequence stars. In Sect. 5, we 
compare our results with observations of supersoft X-ray 
sources and derive clues which may help to identify the 
remaining main sequence star in a supernova la remnant. 
Our conclusions are given in Sect. 6. 



2. Computational method and physical 
assumptions 

The numerical models presented in this work are com- 
puted with a binary stellar evolution code developed by 
Braun (1997) on the basis of a single star code (Langer 
1998, and references therein). It is a l-dimensional implicit 
Lagrangian code which solves the hydrodynamic form of 
the stellar structure and evolution equations (Kippenhahn 
& Weigert 1990). The evolution of two stars and, in case 
of mass transfer, the evolution of the mass transfer rate 
and of the orbital separation are computed simultaneously 
through an implicit coupling scheme (see also Wellstein 
& Langer 1999, Wellstein et al. 1999), using the Roche- 
approximation in the formulation of Eggleton (1983). To 
compute the mass transfer rate, we use the description 
of Ritter (1988). The stellar models are computed using 
OPAL opacities (Iglesias & Rogers 1996) and extended 
nuclear networks including the pp I, II, and HI chains, 
the four CNO-cycles, and the NeNa- and MgAl-cycles (cf. 
Arnould et al. 1999). 



2.1. The mass accretion rate of the white dwarf 

In order to compute the evolution of close main sequence 
star-white dwarf pairs, we invoke the following assump- 
tions. In the computation of the binary system the white 
dwarf is approximated as a point mass (however, cf. 
Sect. 2.3), while the main sequence star is resolved with 
typically 1000 grid points. The systems are started at 
time t:=0 with a zero age main sequence star of mass 
Mi.i and an arbitrary white dwarf mass A^fwD i 
arbitrary orbital separation dj. For most models, we use 
-^WD.i = 0.8 Mq or 1.0 Mq. For dj we consider only val- 
ues which lead to mass transfer during the core hydrogen 
burning phase of the main sequence star, i.e. so called 
Case A mass transfer. 

The mass of the white dwarf (the point mass) is al- 
lowed to vary in accordance with critical mass trans- 
fer rates which were taken from the literature as fol- 
lows. For mass transfer rates M > Mh('^^wd) and M > 
MueiMwo), we allow the white dwarf mass to increase. 
Here, Mu and Mhc are the critical accretion rates above 
which H- or He-burning proceeds such that violent nova 
flashes and consequent mass ejection from the white dwarf 
are avoided. We adopt Mne = 10~^ Mq yr~^ for models 
with a metalhcity of Z = 2%, and 4 IQ-^ Mq yr'^ for 
Z = 0.001 (Fujimoto 1982, Nomoto & Kondo 1991). For 
Mh we rely on Figure 5 of Kahabka & van den Heuvel 
(1997). For Mhc > M > Mh we assume the white dwarf 
mass to grow as well, but by accumulating a degener- 
ate thick helium layer. For M < MEdd := -^Edd/^ and 
M < Mrg, wc assume Mwd = M. Here, Mrg is the crit- 
ical accretion rate above which the white dwarf is assumed 
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to expand to red giant dimensions (Nomoto & Kondo 
1991), and 



L 



Edd 



0.2* (1 + X) 



(1) 



is the Eddington luminosity of the white dwarf, using 
0.2*(1 + X) as the opacity coefficient due to electron 
scattering with X being the hydrogen mass fraction. The 
quantity e = 7 10^^ erg g~^ gives the approximate amount 
of energy obtained per gram of hydrogen burnt into he- 
lium or carbon/oxygen. For larger mass transfer rates we 
assume that the white dwarf has a wind which carries the 
excess mass away (Hachisu et al. 1996). We stop our cal- 
culations for models with Af^ind > 3A/Edd- 

Our assumptions concerning the critical accretion rates 
are similar to those of Li & van den Heuvel (1997). How- 
ever, we deviate from them by adopting a maximum 
possible wind mass loss rate of Afwind > 3AfEdd- For 
M = 2A/Edd, the wind momentum AfwindWoo = A/Edd^oo 
is of the order of the photon momentum LEdd/c. More 
specifically, Afwind = A^Edd implies 



A/wind ^oo '^ocC 

L/c e 
which, for 



2GA/. 



WD 



RwD 



(2) 



(3) 



is of order unity. Our restriction implies that the winds 
we invoke remain in a regime where the wind efficiency 
is undisputed (cf. Lamers & Cassinelli 1999). It lim- 
its our mass loss to rates well below those allowed 
by Li & van den Heuvel (1997). I.e., with AfEdd — 
3.3 10~^(AfwD/ Mq) M0 yr~^, our upper limit is of the 
order of 10"^ M© yr'^ rather than 10"'' Mq yr'^ 

In this context we note that there is also an energy 
limit to radiation driven winds such that M < A/eii ■— 
L/vl^; i.e. 



AfEn = 



L 



-^Edd 



-'escape 



27rci?wD 
0.2*(1 + X)' 



(4) 



This assumes spherical symmetry and ignores the thermal 
energy of the wind. It implies that a star with M = A'/eh 
is invisible, as all the photon energy is used to drive the 

, this results in 



"^escape 



wind. With L = L^dd and Wqo 
AfE„ = 61O-6(i?wD/O.OlR0) Mq yr-i. 

Kato & Iben (1992) and Kato & Hachisu (1994) have 
worked out a theory for optically thick winds which al- 
lows to obtain mass loss rates which can carry well above 
100 times the photon momentum. This has been used by 
Hachisu et al. (1996) and Li & van den Heuvel (1997). 
However, such optically thick winds from stars near the 
Eddington limit may involve processes which limit the 
wind efficiency. For example, if the underlying star car- 
ries noticeable amounts of angular momentum it may 



reach critical rotation before reaching the Eddington limit 
(Langer 1997, 1998). This means that only a fraction of 
the stellar surface, around the equator, is experiencing the 
critical outflow condition, rather than the whole stellar 
surface. The reality of this phenomenon is demonstrated 
by the highly non-spherical, axially symmetric nebulae 
around Luminous Blue Variables (Nota et al. 1995), whose 
outbursts are likely driven by super-Eddington winds (cf. 
Langer 2000). Also nova winds and outflows from super- 
soft X-ray sources are known to be highly anisotropic. 
Another limiting factor is convection or turbulent energy 
transport, which can be very important in Eddington flows 
(Heger & Langer 1996, Owocki & Gayley 1997, Langer 
1997). Furthermore, authors who compute the driving 
force in winds of Wolf-Rayet stars using detailed non- 
equilibrium atomic physics in order to compute the pho- 
ton trapping in the optically thick parts of the wind flow, 
rather than relying on the continuum approximation |^ 
flnd that only about 5% of the stellar photon luminos- 
ity is converted into kinetic wind energy (Lucy & Ab- 
bott 1993). Were this number valid for white dwarfs at 
the Eddington limit, it would imply a mass loss rate of 
M = 0.05AfE = 31O-^(i?wD/O.OlR0) M© yr"!. 

The consequences of our conservative assumption on 
the white dwarf wind efficiency, which deviates from as- 
sumptions on wind mass loss rates in previous studies (cf. 
Hachisu et al. 1996, Kobayashi et al. 1998), are discussed 
in Sect. 3.2. Self-excited wind as proposed by King & van 
Teeseling (1998) are not considered here. 

Other than Li & van den Heuvel (1997), we do not 
allow for the possibility of the partial mass ejection in 
case of weak shell flashes, as we think that the cur- 
rent uncertainties (cf. Prialnik & Kovetz 1995) may not 
make such sophistication worthwhile but rather compli- 
cate the understanding of the obtained results. Our value 
of A/h = 10"*^ Mq yr^i for Z = 2% is in agreement 
with the general conclusion of Prialnik & Kovetz that 
M ^ 10~® Mq yr~^ leads to growing white dwarf masses. 
We also do not include a reduction of the mass accumula- 
tion efficiency of the white dwarf due to winds excited by 
helium shell flashes (Kato & Hachisu 1999). However, in 
order to study the influence of the threshold value for mass 
accumulation on the white dwarf, we investigate the effect 
of an increase of this value by one order of magnitude, as 
outlined in Sect. 3.2. 

2.2. Further white dwarf properties 

Basic properties of the white dwarf are estimated as fol- 
lows: Its radius i?wD is given by the mass-radius relation 



(5) 



^ We note in passing that in order to derive the driving force 
from the continuum approximation, the flux-mean opacity co- 
efflcient needs to be used, not the Rosseland mean opacity, 
which is often used instead. 



4 



N. Langer et al.: Main sequence star + white dwarf binaries as SN la progenitors 



with / = 

9.03101^ em 0-1/3 



2 



cmg"'" (Kippcnhahn & Weigert 1990; see 
also Nauenberg 1972, Provencal et al. 1998). 

We assume that the radiation from the white dwarf 
can be approximated by a Black Body with an effective 
temperature 



'■WD 



(6) 

cm~^ being the Stefan- 



with a = 5.67 10"^ erg s" 
Boltzmann constant. 

We estimate the white dwarf's angular momentum and 
rotation frequency by assuming it to be zero initially, and 
that the specific angular momentum of the accreted mat- 
ter is j = WKepler-RwD, with ^Kepler = -\/GMwd/^?WD- 

The equation for the accumulated angular momentum 



J:= f Mwujdt' 
Jt'=o 

yields 



4/3 
WD 



4/3 
WD,i 



(7) 



(8) 



Assuming then rigid rotation for the white dwarf interior 
(cf. Kippenhahn 1974) we can estimate its angular velocity 
as 

J 



(9) 



where k is the dimensionless radius of gyration . The ratio 
O = w/wKepier bccomes 




/ MwD,i 
V ^WD 



4/3N 



(10) 



(cf. Papaloizou & Pringle 1978, Livio & Pringle 1998). The 
implications of this relation are discussed in Sect. 4.1. 

2.3. Single star evolution 

We have constructed stellar model sequences for two 
metalhcities, Z = 0.02 and Z = 0.001. The initial helium 
mass fraction is computed as F = Ypm + (dY/dZ)Z, us- 
ing Ypni = 0.24 as primordial helium mass fraction and 
dY/dZ = 2. The resulting values are Y = 0.280 and 
Y = 0.242 for the high and low mctallicity considered 
here, respectively. The relative abundances of the met- 
als are chosen according to the solar system abundances 
(Grcvcssc & Sauval 1998). We have computed all mod- 
els with extended convective cores ("overshooting") by 
0.2 pressure scale heights (aover = 0.2). The resulting 
tracks in the HR diagram arc very similar to those of 
Schaller et al. (1992), that of our Z = 0.02 models for 
2Mq and 1.7 M© and of our 2M0 sequence at Z = 0.001 

— the only three sequences which we can directly compare 

— are virtually identical. 




Fig. 1. Evolution of the stellar radius as a function of time 
for single stars of about solar mctallicity (Z=0.02) in the 
mass range from 1.7 to 2.4 computed with convective 
core overshooting (aovcr = 0.2), from the zero age main 
sequence until shortly after core hydrogen exhaustion. 

Table 1. Comparison of surface properties and of the 
Kelvin-Helmholtz time scale tkh.tams = GM"^ /{2RL) for 
models of 2 Mq sequences computed with metalhcities of 
Z = 0.02 and Z = 0.001, at the zero age main sequence 
(ZAMS) and the terminal age main sequence (TAMS) 





Z = 0.02 


Z = 0.001 


izAMS 


15.9 Lo 


26.3 Lq 


Toff, ZAMS 


9158 K 


12474 K 


JJzAMS 


1.60 Rq 


1.11 Ro 


TKH,ZAMS 


2.5Myr 


2.2Myr 


-I/TAMS 


22.3 Lc3 


56.1 L0 


Teff.TAMS 


6600 K 


10452 K 


-Rtams 


3.65 R0 


2.31 R© 


Tkh.tams 


0.75Myr 


0.48Myr 



Figs. 1 and 2 show the time evolution of the radii dur- 
ing core hydrogen burning for our single star models at 
the two metalhcities considered in this work, for the mass 
range which is relevant in the context of this paper. Dur- 
ing core hydrogen burning, the radii increase by factors 
2.0... 2. 5 and 1.7... 2. 3 for the higher and lower considered 
metallicity, respectively, with larger values corresponding 
to larger masses. 

The radii of the metal poor stars are nearly a factor 
of 2 smaller than those of stars with a comparable mass 
and evolutionary stage at Z = 0.02. This has consequences 
for the binary evolution models discussed below, i.e., the 
orbital periods in Case A systems are much smaller for 
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Fig. 2. Same as Figure 1, but for fow metallicity stars 
(Z=0.001) in the mass range 1.4 to 2Mq. 

smaller metallicity. Furthermore, the metal poor main se- 
quence stars are hotter and more luminous compared to 
stars of the same mass a.t Z = 0.02. This is in agreement 
with previous models of stars of comparable masses and 
metallicities (e.g., Schaller et al. 1992). Table 1 gives the 
quantitative details of models from our 2 sequences at 
Z ~ 0.02 and Z ~ 0.001 at the beginning and at the end 
of core hydrogen burning. Although the models at lower 
metallicities are more compact, they are also more lumi- 
nous than the metal richer models, with the consequence 
of a shorter Kelvin-Helmholtz time scale. This has con- 
sequences for the mass transfer rates, as discussed in the 
next section. 

2.4- Binary evolution: examples 

To illustrate our approach, we consider our System No. 6. 
Initially, it consists of a 2 Mq zero age main sequence star 
and a 1 Mq white dwarf (treated as point mass) in a cir- 
cular orbit with a separation of di = 7.19 R©. According 
to Kepler's third law, the initial period Pi is 



' VGV MMS,i + MwD,i 



(11) 



i.e.. Pi = 1.29 d in this case. With this initial set-up, we 

skip the previous evolution of the system, i.e. the evolution 
of the white dwarf progenitor and the common envelope 
and spiral-in phase which brought the two stars close to- 
gether. 

We set the time f = at core hydrogen ignition of our 
main sequence star. In principle, this neglects the dura- 
tion of the previous evolution of the system, i.e. the evo- 
lutionary time of the white dwarf progenitor. However, as 
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Fig. 3. Evolution of the orbital period as function of time 
for System No. 6 (cf. Table 1). The curve starts at the 
onset of mass transfer at t ~ 1.166 Gyr. At t ~ 1.169, 
the white dwarf mass has grown to 1.4 M©; this time is 
marked by an asterisk on the curve. The dashed part of 
the curve shows the continuation of the orbital period evo- 
lution assuming that the white dwarf does not perform a 
supernova explosion. 

we shall see below we deal with rather high initial white 
dwarf masses, i.e. relatively massive white dwarf progen- 
itor stars (^ 5M0). Therefore, as the white dwarf starts 
to accrete at a system ago of the order of the evolutionary 
time scale of the donor star (1.4...2.4 M©; see below) the 
so defined time yields a good estimate for the age of the 
systems at the time of the supernova explosion (cf . Umeda 
et al. 1999). 

In our System No. 6, mass transfer starts at t 
1.166 Gyr, at a central hydrogen mass fraction of the main 
sequence star of Xc — 0.17. The radius of the star has then 
grown from Pi ~ 1.60 R© at the ZAMS to P ~ 3.16 Rq, 
as wc! use Eggleton's (1983) approximation for the Roche 
radius of the main sequence star 



Pl 



0.49(7' 



2/3 



0.6g2/3 ln(l-hgi/3) 



(12) 



with q := Mms/-^wd- 

During the initial phase of the mass transfer evolu- 
tion, the mass transfer rate M is still smaller than Mh (cf. 
Sect. 2.1), and the white dwarf mass can not increase. In- 
stead, all accreted mass is assumed to be lost in nova out- 
bursts, carrying away the specific orbital angular momen- 
tum of the white dwarf, leading to a decrease of the orbital 
separation (cf. Podsiadlowski et al. 1992). The amount of 
mass lost during this phase, and consequently the change 
of the orbital parameters, is quite insignificant in most 
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1.174 



Fig. 4. Evolution of tlie mass transfer rate as function of 
time for System No. 6 (cf. Table 1). The curve starts at 
the onset of mass transfer at t ~ 1.166 Gyr. At t ~ 1.169, 
the white dwarf mass has grown to 1.4 M©; this time is 
marked by an asterisk on the curve. The dashed part of 
the curve shows the continuation of the mass transfer rate 
evolution assuming that the white dwarf does not perform 
a supernova explosion. 

cases. E.g., in our example 0.003 are transferred and 
lost in nova outbursts. However, we emphasise that the 
time scale of this first phase may be non- negligible, as it 
may be of the same order of magnitude as the major mass 
accretion phase (cf. Sect. 3). 

When the mass transfer rate exceeds the critical rates 
for hydrogen and helium burning (Sect. 2.1), we allow the 
white dwarf mass to grow. As angular momentum con- 
servation again leads to a shrinkage of the orbit as long 
as A^MS > Af\vD: the mass transfer is thermally unsta- 
ble. Since the mass-radius exponents of our main sequence 
stars are positive, i.e., mass loss leads to smaller radii (Rit- 
ter 1996), the mass transfer is stabilised due to the thermal 
disequilibrium of the main sequence star, and the resulting 
mass transfer rates are of the order of magnitude of 



M ~ (MMS,i - MwD,i)/T-KH 



(13) 



(e.g., Rappaport et al. 1994; however, see Sect. 3.1). 

Fig. 3 shows the evolution of the orbital separation 
with time for System No. 6 throughout the phase of 
the thermally unstable mass transfer. Within our ap- 
proximations, the white dwarf mass reaches 1.4 Mq at 
t ~ 1.169 Gyr. Although it is likely that the white dwarf 
would explode roughly at this point (see below), we fol- 
low the further evolution of the system ignoring this, for 
several reasons. Most important, the continued evolution 
allows us to estimate how changes in our basic assump- 
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Fig. 5. Evolution of the mass transfer rate as function of 
time for System No. 61 (cf. Fig. 8). 
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Fig. 6. Evolution of the orbital separation as function of 
time for System No. 61. 

tions might affect the fate of the white dwarf. E.g., since it 
appears to be undisputed that the white dwarf mass can 
grow for mass transfer rates M ^ 10~^ M© yr~^ (Nomoto 
& Kondo 1991, Prialnik & Kovetz 1995, Kahabka & van 
den Heuvel 1997), we define Mr := / Mrdt with 



M7 = 



for M < IQ-"^ Mq yr-i 
MwD for M > 10"'^ Mq yi'^ 



(14) 



I.e., My gives the amount of mass by which the white 

dwarf grows at accretion rates above 10~^ Mq yi~^, in- 
cluding the continued evolution beyond Mwd = 1-4 Mq. 
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For System No. 6, the point mass has grown to 1.75 Mq 
by the time the mass transfer rate falls below the critical 
rate (~ 10~* Mq yr~-^). The quantity My in this system is 
M7 = 0.60 M0, which means that were the critical accre- 
tion rate as high as 10~^ Mq yr^^, the point mass would 
still have grown to 1.60 M©. This can be understood from 
Fig. 4, which shows the mass transfer rate as function of 
time for System No. 6. It can be seen that the mass trans- 
fer rate remains above 10~^ Mq yr~^ for several million 
years after the potential supernova explosion. 

Note that this estimate is not fully self-consistent. I.e., 
were all the accreted mass lost through nova outbursts as 
long as the mass transfer rate were below 10-"^ Mq yr , 
the orbital evolution would differ from that of our model. 
However, as the mass loss from the white dwarf would keep 
the mass ratio q := MmsM^wd above one for a longer 
time, the orbit would keep shrinking for a longer time, 
which would keep the mass transfer rate higher than in 
our model (cf. Sect. 3.1). Thus, the value of M7 which 
we derive is in fact a lower limit to the mass which is 
transfered at rates above 10"'' Mq yr~^. 

Comparing Figs. 3 and 4 shows, that the mass transfer 
rate drops to small values only several million years after 
the minimum orbital separation is achieved. The reason 
is that at the time of minimum separation the main se- 
quence star is still more compact than its thermal equi- 
librium configuration. I.e., even though the orbit does not 
shrink any more, the main sequence star expands towards 
its thermal equilibrium radius and drives further mass 
transfer thereby. 

In Figs. 5 and 6 we show the long-term evolution of the 
mass transfer rate and of the orbital period for another 
system, No. 61, initially consisting of a low metallicity 
1.8 Mq main sequence star and a 1 Mq white dwarf or- 
biting with a period of 0.70 d. One can see the thermally 
unstable mass transfer phase lasting for several million 
years, beyond which mass transfer continues only on the 
nuclear time scale of the main sequence star, i.e., several 
10® yr. Consequently, the mass transfer rate drops to some 
10"^'^ Mq yr"'^. The system then resembles a Cataclysmic 
Variable, evolving on a time scale of 10®...10^'' yr. Note 
that in CVs, this time scale may become shorter due to 
angular momentum loss through magnetic braking (Vcr- 
bunt & Zwaan 1981) which is ignored in the present study. 
Magnetic braking is not relevant for the supernova la pro- 
genitor evolution for two reasons. First, the time scale of 
the thermally unstable mass transfer is only of the or- 
der of several million years, which is too short to allow 
a significant amount of angular momentum loss through 
this mechanism. Second, our main sequence stars do de- 
velop convective envelopes only in the final phase of the 
thermally unstable mass transfer phase. I.e., most of the 
time they have radiative envelopes and thus supposedly 
no magnetic wind. For the study of the long term evolu- 
tion of those systems which fail to bring the white dwarf 
to explode as a supernova, magnetic braking might be rel- 



evant. This is, however, beyond the scope of the present 
investigation. 

3. Mass treinsfer and white dwcirf evolution 

3.1. Mass transfer rates 
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Fig. 7. Evolution of the mass transfer rate as function 

of time for systems No. 0, 2, 3, 7 and 13, which have 
a metallicity of Z=0.02 and white dwarf companions with 
initially 1 Mq (see also Fig. 8). System No.O (leftmost line) 
is stopped when the mass transfer rate exceeds the allowed 
upper limit for the wind mass loss rate (cf. Sect. 2.1). 
For the other four curves, higher peak mass transfer rates 
correspond to larger initial main sequence star masses 
(2.3 Mq, 2.1 Mq, 2.0 Mq, and 1.8 Mq). The time t = 
is defined by the onset of mass transfer. The star sym- 
bol indicates the time when the white dwarf has reached 
1.4 Mq. Beyond that point, the graphs are continued as 
dotted lines. 

In this Section, we deal with the mass transfer rates M , 
and we emphasise that the mass accumulation rate of the 
white dwarf Mwd may be smaller than the mass transfer 
rate if the latter is above or below the threshold values 
defined in Sect. 2.1. We want investigate the dependence of 
the mass transfer rate, and of its time dependence, on the 
various initial parameters of our binary systems. Although 
Eq. (13) gives the order of magnitude of the mass transfer 
rate during the thermally unstable phase, we will see that 
it fails to reproduce all the physical dependences correctly. 

First consider the dependence of the mass transfer rate 
on the initial mass of the main sequence component Mms.i- 
Figures 7 and 8 show the mass transfer rate as function 
of time for systems with white dwarf initial masses of 
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Table 2. Key properties of interacting main sequence star + white dwarf systems with Z = 0.02. The cohimns have 
the following meanings. (1) system number, (2) main sequence star initial mass, (3) white dwarf initial mass, (4) initial 
orbital period, (5) minimum period, (6) maximum possible CO-mass (see text), (7) main sequence star mass when 
■^WD = 1.4 M0, (8) total mass loss due to winds (9) see Eq. 13, (10) maximum mass transfer rate (11) maximum 
X-ray luminosity of the white dwarf (12) core hydrogen mass fraction of main sequence star at onset of mass transfer 
(13) system age when Mwd = 1-4 Mq 
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-^WD,i = 1 Mq and various initial main sequence masses, 
for Z = 0.02 and Z = 0.001, respectively. For both metal- 
licities, there is a clear trend to larger maximum mass 
transfer rates for more massive main sequence stars caused 
by the shorter thermal time scale of more massive main 
sequence stars. Eq. (13), with tkh := GM^ /{2RL), does 
reproduce the maximum mass transfer rates within 30% 
for all sequences shown in Fig. 8. However, it overesti- 
mates those of the sequences shown in Fig. 7 by factors 
3... 8, larger values corresponding to smaller initial main 
sequence star masses. 

In all our models, there is a time delay from the onset of 
the mass transfer (defined as < = in Figs. 7 and 8) to the 
time when the mass transfer rate has grown sufficiently 
to allow the white dwarf mass to grow. This delay is of 
the order of the thermal time scale of the main sequence 
star, i.e. it is longer for smaller masses. We emphasise 
that our method to compute the mass transfer rate (Ritter 
1988, see also Braun 1997) allows its reliable computation 
also for the beginning and the end of the mass transfer 
evolution. Assuming a nova outburst would occur after the 
accumulation of ~ 10~^ Mq (cf. Prialnik & Kovetz 1995) 
and mass accretion rates of the order of 10~^ Mq yr~^ 
— i.e. nova recurrence times of about 10^ yr — implies 
of the order of thousand nova outbursts in our typical 



supernova la progenitors before the white dwarf mass can 
start growing. 

Fig. 9 shows the mass transfer rate as a function of 
time for three systems with the same initial main sequence 
star mass but with different initial white dwarf masses. 
The time delay from the onset of the mass transfer until 
the mass transfer rate exceeds ^ 10~^ Mq yr~^ and the 
white dwarf mass can start growing is very similar for all 
three systems. The delay is determined by the thermal 
time scale of the main sequence star. 

Most striking in Fig. 9 is the feature that much larger 
mass transfer rates are achieved for smaller white dwarf 
masses. Although this trend is also expected from Eq. (13), 
the order of magnitude of the effect seen in Fig. 9 is much 
larger than what Eq. (13) predicts. We find that, although 
initially less massive white dwarfs need to accrete more 
mass to reach the Chandrasekhar mass, the best super- 
nova la candidate system of those displayed in Fig. 9 may 
actually be the one with the smallest initial white dwarf 
mass. This is so since the same donor star transfers much 
more mass for smaller initial white dwarf masses, and that 
even at higher mass transfer rates. 

For otherwise fixed system parameters, more mass is 
transfered for smaller initial white dwarf masses since the 
minimum orbital separation — coincides with the time 
when Mms = .^wd in the conservative case — is ob- 
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Fig. 8. Evolution of tlie mass transfer rate as function 
of time for low metallicity systems No. 58, 61, 63, 65, 
67 and 70. The initial white dwarf mass is 1 M0 for all 
six cases. Higher peak mass transfer rates correspond to 
larger initial main sequence star masses, from 1.9 to 
1.4 Mq in steps of 0.1 Mq. System No. 58 (leftmost line) is 
stopped when the mass transfer rate exceeds the allowed 
upper limit for the wind mass loss rate (cf. Sect. 2.1). 
Note that the scale of the vertical axis is identical to that 
in Fig. 7, but the represented initial masses of the main 
sequence components is lower. 



tained only after more mass is transferred. Higher mass 
transfer rates arc achieved since smaller minimum orbital 
separations are obtained for smaller values of MwD,i- For 
conservative evolution of a given binary system, the or- 
bital separation d can be expressed as 



d=J 



^Mms + Mwd 



where 

^ 27rd2 Mms^wd 



P M, 



MS 



WD 



(15) 



(16) 



is the constant orbital angular momentum. Therefore, a 
given initial separation di relates to the minimum orbital 
separation dmin as 



d„ 



di 



= 4 



MMS,iMwD,i 
'(MMS,i + MwD,i) 



(17) 



and as for fixed system mass the period and separations 
are related as oc d'^ it is 



Pi 



= 4- 



.^MS,i-^WD,i 
" (MMS.i + MwD,i)^ 



(18) 



Fig. 9. Evolution of the mass transfer rate as function 
of time for systems No. 13, 39 and 53, i.e. systems with 
a metallicity of Z = 0.02 and with initial main sequence 
star masses of 1.8 Mq, but with three different initial white 
dwarf masses, as indicated. (Cf. also caption to Fig. 7.) 



I.e., let us consider two main sequence stars of the same 

mass, starting to transfer mass onto their white dwarf 
companions at the same orbital separation di. The min- 
imum separation will be smaller for the binary with the 
smaller initial white dwarf mass, say System A. Since the 
radius of main sequence stars in the considered mass range 
decrease for increasing mass loss rates, the mass loss rate 
of the main sequence star in System A — i.e., its mass 
transfer rate — needs to be larger in order to fit the main 
sequence star into a smaller volume. 

The fact that the white dwarf in System No. 53, which 
has an initial mass of 0.7 M©, does not reach 1.4 M© but 
only 1.37 M© is due to the fact that during the peak of 
the mass transfer the rate slightly exceeds the Eddington 
accretion rate MEdd (which is smaller for smaller white 
dwarf masses; cf. Sect. 2.1), and this system loses 0.27 M© 
to a wind. The other two systems shown in Fig. 9, No. 13 
(MwD,i = IM0) and No. 39 (MwD,i = 0.8 Mg), which 
avoid winds, can grow the CO-white dwarf to 1.45 Mq and 
1.46 Mq, respectively. In fact, the system with the largest 
initial white dwarf mass. System No. 13, is least likely to 
produce a Type la supernova, since in this system all mass 
is transferred at rates below 10~^ Mq yr"'^. I.e., My = 
in this case, while My = 0.45 Mq for System No. 39, and 
My = 0.51 Mq for System No. 53 (cf. Table 2). 

I.e., the effect that systems with smaller initial white 
dwarf masses are better Type la supernova progenitor 
candidates is only limited by the smaller upper limits to 
the white dwarf mass accumulation rate for smaller white 
dwarf masses. 
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Wc want to point out that, for the Case A systems con- 
sidered in this work, the maximum mass transfer rate can 
vary by up to a factor of 2 as function of the initial period 
(cf. Tables 2 and 3). One would expect larger mass trans- 
fer rates for initially wider systems, since in this case the 
main sequence star is more extended and more luminous 
at the onset of the mass transfer, and thus has a shorter 
Kelvin-Helmholtz time scales (cf. Table 1). This expecta- 
tion, which is also reflected in Eq. (13), is fulfilled rather 
well for most of our low metallicity systems (cf. Table 3) . 
However, at Z = 0.02 we find mostly decreasing maximum 
mass transfer rates for increasing initial periods and oth- 
erwise fixed initial system parameters (cf. Table 2). This 
means that Eq. (13) can not be used to predict trends of 
the mass transfer rate as function of the initial period or 
the system metallicity, and shows the limitations of sim- 
plifying approaches to the study of accreting white dwarfs 
in binary systems. 



Fig. 10. Evolution of the mass transfer rate as function of 3.2. Evolution of the white dwarf mass 
time for systems No. 9 and 60, which have identical initial 

main sequence star and white dwarf masses, but different 

metallicities, as indicated. (Cf. also caption to Fig. 7.) 
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The dependence of the mass transfer rate on the metal- 
licity of the main sequence star is elucidated in Fig. 10. It 
shows two systems with identical initial main sequence 
star and initial white dwarf masses but with different 
metallicities. Wind mass loss is negligible in both cases. 
The difference in the maximum mass transfer rate of Sys- 
tems No. 9 and No. 60 — which both start with a main 
sequence star of 1.8 Mq and a white dwarf of 1 Mq — is 
a factor of four. Both systems have initial periods close 
to the shortest possible initial period. This large differ- 
ence is not due to different stellar Kelvin-Helmholtz time 
scales tkh- Although the low metallicity stars are more 
luminous, they are also more compact both effects on tkh 
almost cancel out (cf. also Table 1). At the onset of mass 
transfer, it is tkh = 2.6 Myr for the main sequence star 
in System No. 9, while the corresponding value for Sys- 
tem No. 60 is Tkh = 2.4 Myr. This is also reflected in the 
similarity of the turn-on times for the mass transfer (cf. 
Fig. 10). 

Low metallicity systems have larger mass transfer rates 
compared to systems with solar abundances (Figs. 7 
and 8). While the range of mass transfer rates covered 
in both figures is the same. Fig. 7 {Z = 0.02) shows 
systems with initial main sequence masses in the range 
2.4...1.8Mq, while those in Fig. 8 {Z = 0.001) are in the 
range 1.9...1.4M0. Based on our detailed models, Eq. (13) 
is valid for low metallicity within a factor of two. For so- 
lar metallicities, the mass transfer rates are systematically 
lower by a factor of 5. I.e., on average the low metallicity 
systems have, for the same initial stellar masses, five times 
higher mass transfer rates than the systems a,t Z = 0.02. 




t/Myr 

Fig. 11. Evolution of the white dwarf mass as function of 
time for systems No. 0, 2, 3, 7, and 13 (cf. Fig. 7.) 

The dependence of the mass transfer rate on various 
parameters discussed in Sect. 3.1 has important implica- 
tions for the evolution of the white dwarfs. Figs. 11 and 12 
illustrate the time evolution of the white dwarf masses 
for the same systems for which the evolution of the mass 
transfer rate M has been displayed in Figs. 7 and 8. We re- 
call that |M| 7^ |Mwd| due to the restrictions on the mass 
accretion rate outlined in Sect. 2. I.e., the white dwarf 
mass can start to grow only 0.5... 3 Myr after the onset of 
the mass transfer, due to the occurrence of nova outbursts 
(cf. also Sect. 3.1). 
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Fig. 12. Evolution of the white dwarf mass as function of 
time for systems No. 58, 61, 63, 65, 67 and 70 (of. Fig. 8.) 
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Fig. 13. Maximum achievable CO-mass Mco (see text; 
solid line and solid dots), Mco plus the total mass lost to 
a wind (dotted line and dots), and initial white dwarf mass 
(IMq for all systems shown here) plus M7 (cf. Eq. (14)) 

— i.e. the at least achieved CO-mass in the white dwarf 

— for the low mctallicity systems No. 58, 61, 63, 65, 67 
and 70 (cf. Figs. 7 and 12, and also Figs. 16 and 17). 



To demonstrate the effect of the upper and lower crit- 
ical accretion rates for the achievable white dwarf masses, 
we have plotted in Fig. 13 the maximum possible CO-mass 
in the white dwarf — ignoring the possible occurrence of 
a supernova event at Mwd — 1.4 Mq — as function of the 



initial main sequence mass for low mctallicity systems with 
an initial white dwarf mass of 1 Mq. The sharp drop of the 
curve at Mms.i = 1-9 M© is due to the fact that the mass 
accretion rate exceeds three times the Eddington accre- 
tion rate of the white dwarf shortly after the onset of the 
mass transfer in the system with Mms.i = 1.9 Mq, which 
we use as criterion to stop the calculations (cf. Sect. 2), 
assuming that the white dwarf would form an extended 
hydrogen-rich envelope and the two stars in the system 
would merge. 
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Fig. 14. Maximum achievable CO-core masses as function 

of the initial mass of the main sequence star, for various 
initial white dwarf masses and for the two metallicities 
considered here, as indicated (cf. also Figs. 13 and 17). 

Fig. 13 also shows the sum of Mco and the total 
amount of mass lost from the system due to a white dwarf 
wind (see also Tables 2 and 3). It indicates that winds, 
and thus the upper critical accretion rate, are unimpor- 
tant for initial main sequence masses Mms.i ^ 1.6 M©, 
at Z = 0.001. As discussed in Sect. 2, the wind effi- 
ciency is quite uncertain, particularly at low mctallicity. 
I.e., Hachisu et al. (1996) and Kobayashi et al. (1998) as- 
sume a much higher wind efficiency compared to our as- 
sumption (Sect. 2), but assume that the white dwarf winds 
break down all together at low mctallicity. A lower wind 
mass loss rate might lead to a merging of our systems with 
Mmsa > 1.6 Mq in Fig. 13, rather than to a Type la su- 
pernova. A higher wind mass loss rate, on the other hand, 
might allow also white dwarfs in systems with donor star 
masses larger than 1.8 Mq to reach the Chandrasekhar 
mass. We note that, as the maximum mass transfer rate 
rises very sharply with increasing initial main sequence 
mass (Figs. 7 and 8), the maximum main sequence mass 
is not very sensitive to the assumptions on the wind mass 
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Fig. 15. CO-core masses achieved if only mass accreted 
with rates above 10"'' M© yr~^ are considered, as function 
of the initial mass of the main sequence star, for various 
initial white dwarf masses and for the two metallicities 
considered here, as indicated (cf. also Figs. 13 and 16). 

loss rate. E.g., Li & van den Heuvel (1997), following 

Hachisu et al. (1996), allowed wind mass loss rates of up 
to 10~^ M0 yr~^, i.e. roughly 100 times more than in 
our calculations. This shifts the maximum main sequence 
mass from 2.3 in our case to 2.6 M0 in their case, for 
Case A systems with an initial white dwarf mass of 1 Mq 
and a metallicity of 2%. 

In order to estimate the effect of a considerably re- 
duced mass accumulation efficiency of the white dwarf due 
to weak hydrogen shell flashes (cf. Prialnik & Kovetz 1995) 
or winds excited by helium shell flashes (Kato & Hachisu 
1999), the lower curve in Fig. 13 shows MwD,i + Mr, i.e. 
the maximum achievable CO mass assuming that mass 
accumulation on the white dwarf is only possible for 
M > 10"'^ M0 yr-i (cf. Eq. (14)). It shows that under 
these assumptions the white dwarfs in the low metallicity 
systems with initial main sequence masses of I.5...I.8M0 
would still be able to grow to 1.4 M0, although not to sig- 
nificantly larger values. Furthermore, we note from Fig. 13 
(cf. also Figs. 14 and 15) that a reduction of the limiting 
mass accretion rate for mass accumulation by one order 
of magnitude has little effect on the upper limit of Mms.i 
in supernova progenitor systems. 

Figs. 14 and 15 show the complete picture of the out- 
come of our study for the achievable white dwarf mass 
as function of the system parameters. While Fig. 14 gives 
the optimistic view, i.e. applying the lower critical accre- 
tion rates defined in Sect. 2, Fig. 15 shows the graphs 
for MwD.i + Afy rather than for Mco- When consider- 
ing Figs. 14 and 15, it is important to keep in mind that 



the variation of the initial period of the considered sys- 
tems would convert each line in these figures into a band 
with an average width of the order of 0.1 Mq (cf. Tables 2 
and 3). These figures allow the following conclusions. 

1. The initial main sequence masses from which Type la 
supernovae at low metallicity can be drawn (~ 
1.45 Mq. .. 1.85 Mq) are much smaller than at Z = 0.02 
(~ I.7M0...2.35M0). 

2. The initial masses of the white dwarfs required for a 
Type la supernova are — on average — about 0.2 Mq 
larger at Z = 0.001 compared to the high metallicity 
systems. I.e., note that in Figs. 14 and 15 the curves 
for MwD,i = 1 Mq at Z = 0.001 are similar to those 
for MwD,i = 0.8 Mq at Z = 0.02, only shifted to lower 
initial main sequence star masses. 

3. Since smaller initial mass ratios MwD,i/-^MS,i lead to 
larger mass transfer rates, more; mass is transfcrcd in 
systems with smaller initial white dwarf masses than 
for lower initial white dwarf masses. This effect may 
giv(^ systems with small initial white dwarf masses 
i.e. perhaps as low as MwD,i = 0.7 Mq — the possibil- 
ity to produce a Type la supernova. 

4. Evolution of further system parameters 
4-1- White dwarf spin 

In Sect. 2.2, we showed that, within simple approxima- 
tions, the spin of the accreting white dwarfs at any given 
time depends only on the amount of matter accreted up 
to that time. From Eq. (10) follows that the largest spin 
at the time of the supernova explosion is expected in those 
systems that start out with the least massive white dwarfs. 
For initial white dwarf masses of 1.2 Mq, 1.0 Mq, 0.8 Mq, 
and 0.7 Mq, the assumption of homogeneous white dwarfs 
(i.e., fc^ = |) leads to ratios of rotational to critical rota- 
tional velocity fl of 0.35, 0.67, 0.98, and 1.13, respectively. 
A more realistic value of fc = 0.4 (Ritter 1985) results 
even in Q = 0.85, 1.68, 2.45, and 2.83. While these num- 
bers should not be taken literally in particular, values 
of n > 1 are of course not plausible — they elucidate 
the possibility that many of the exploding white dwarfs 
in Type la supernovae may be rotating at a speed close 
to break-up. According to Figs. 14 and 15 it is not ex- 
cluded that white dwarfs with initial masses of 0.7 Mq 
may contribute to the Type la supernovae, even at low 
metallicity. 

This point of view is at least partly supported by ob- 
servational evidence. While isolated white dwarfs appear 
to rotate very slowly {viot ^ 50kms~^; Heber et al. 1997, 
Koester et al. 1998) — which is also expected from recent 
single star models with rotation (Langer et al. 1999) — 
those in CVs can be much larger, i.e. up to 1200 kms~^ 
(Sion 1999). As the white dwarfs in CVs are accreting, 
this shows that accreting white dwarfs are in fact spun-up. 
The fact that the white dwarfs in CVs are not spinning as 
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Fig. 16. Evolution of the orbital separation as function of 
time for systems No. 4, 5, 6, 7, and 8. The time t = is 
defined by the onset of mass transfer. A star symbol indi- 
cates the time when the white dwarf has reached 1.4 M©. 
Beyond that point, the graphs are continued as dotted 
lines. 

rapidly as expected from Eq. (10) is interpreted by Livio & 
Pringle (1998) as being due to angular momentum loss in 
nova explosions which must occur in typical CV systems. 

From our simple approach, we expect that by the time 
the white dwarf mass gets close to the Chandrasekhar 
limit, it may rotate with a significant fraction of the break- 
up velocity. First polarisation studies of la supcrnovae 
seem to indicate that the degree of polarisation in the 
supernova spectra is very low (Wang et al. 1996), which 
makes a strongly deformed white dwarf as initial configu- 
ration for the explosion rather unlikely. A confirmation of 
these results on a solid statistical basis would imply that 
cither typical initial white dwarf masses arc rather high, 
or that the white dwarfs can lose angular momentum dur- 
ing their accretion phase through mechanisms yet to be 
found. 

4-2. Orbital evolution 

In all systems the mass transfer rate remains initially on a 
low level for about one Kelvin- Helmholtz time scale of the 
main sequence star. Nova outbursts are to be expected 
for the first 0.5...2Myr of the mass transfer evolution. 
This is simulated as a continuous process in our calcu- 
lations, where all transferred mass leaves the system im- 
mediately, carrying the specific orbital angular momentum 
of the white dwarf. Both, mass transfer and nova winds 
at this stage (Mus > ^wd) lead to a decrease of the 
orbital separation. The effective shrinkage of the orbit is. 



however, small since the amounts of mass involved in the 
mass transfer and winds during this phase are small (cf. 
Sect. 3). We have not included the possibility of frictional 
angular momentum loss during nova outbursts (e.g., Livio 
& Pringle 1998), as this could have been compensated by 
a slight increase of our initial periods which are treated as 
a free parameter anyway. 

Once the mass transfer rate exceeds the critical rates 
Mh and Mue (cf. Sect. 2), we assume the mass transfer to 
be conservative as long as M < MEdd (Eq. (1)). In that 
case, mass and angular momentum conservation lead to a 
shrinkage of the orbit as long as Mms > AfwD, and to a 
widening thereafter. Fig. 16 shows the evolution of the or- 
bital separation with time for five systems with the same 
initial main sequence star and white dwarf mass but with 
different initial periods. As none of these systems devel- 
ops a super-Eddington wind, the minimum orbital sepa- 
rations follow from Eq. (18) (neglecting the nova winds). 
As Eq. (18) can be expressed as 

^^^■^'(^+1?) ' ^^^^ 

using Qi = MwD,i/MMS,i, shorter minimum periods are 
achieved for shorter initial periods, and for smaller initial 
mass ratios qi < 1. 

Those systems which evolve a super-Eddington wind 
can evolve to considerably shorter periods than the con- 
servative systems (Eq. (19)). The main reason is that the 
condition Mms > ^wd remains fulfilled for a much longer 
time than in the conservative case. E.g., Systems No. 1 
and No. 2, which lose about 0.28 and 0.39 to a 
wind, evolve to minimum periods of 6.5 h and 21.6 h, re- 
spectively (cf. Table 1), while their minimum periods in a 
conservative evolution would have been 7.5 h and 25.2 h. 

4-3. Evolution of the main sequence star 

A discussion of the properties of the main sequence stars in 
the supernova la progenitor systems presented before may 
be interesting for two reasons. First, it might be observ- 
able during the accretion phase, where our systems might 
appear as supersoft X-ray sources. Second, the main se- 
quence star may survive the explosion of the white dwarf 
and may then serve as an observable witness of the super- 
nova progenitor evolution. 

To elucidate the first point, we have plotted in Fig. 17 
the evolutionary track of the main sequence component 
of System No. 30 in the HR diagram, in comparison to 
normal single stars tracks for comparable initial masses. 
It is evident that this star, once the Roche lobe overfiow 
sets in, reduces its luminosity significantly, roughly by a 
factor of ^ 30. However, from Fig. 17 it is not clear which 
fraction of the luminosity decrease is due to the fact that 
the stellar mass of the main sequence star becomes smaller 
— from 2.1 Mq to about 0.8 at the time when the 
white dwarf has reached 1.4 Mq — , and which fraction 
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Fig. 17. Evolutionary track in the HR diagram of the 
main sequence component of a 2.1 M0+O.8 M0 main se- 
quence star + white dwarf binary with a metallicity of 2% 
and an initial period of 1.23 d (System No. 30), from zero 
age until the supernova explosion of the white dwarf (thick 
solid line) . Dotted lines correspond to evolutionary tracks 
of single stars of 2.1, 1.9, and 1.7 M©. The two dashed 
lines denote the zero age main sequence (left line) and the 
terminal ago main sequence (corresponding to core hydro- 
gen exhaustion). The star symbol denotes the time when 
the white dwarf reaches 1.4 Mq. 

is due to the deviation of the star from global thermal 
equilibrium during the rapid mass loss phase. 

In order to understand which luminosities the main se- 
quence components of our systems can achieve in general, 
and to what extent the reduction of the luminosity can 
be understood in terms of the mass reduction, we show in 
Figs. 18 and 19 the mass-luminosity evolution for a sam- 
ple of our systems in comparison with the mass-luminosity 
relation for normal main sequence stars, for Z = 0.02 and 
Z = 0.001, respectively. These figures show that in most 
systems the main sequence components are significantly 
underluminous with respect to their mass. 

The inspection of Figs. 18 and 19 reveals the following 
features. First, the systems with the shorter initial peri- 
ods, which start out at the lower boundary of the main 
sequence band in these figures, evolve to the zero age main 
sequence position corresponding to their final mass once 
the mass transfer rates become small and thermal equilib- 
rium is restored. In real binaries the white dwarf might 
explode before this happens, but the models shown in 
Figs. 18 and 19 are followed up to the end of the ther- 
mally unstable mass transfer phase, which corresponds to 
the endpoint of the dotted parts of the tracks. This result 
is easily understood, since when the mass transfer starts 
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Fig. 18. Evolutionary tracks of the main sequence stars 
in the mass-luminosity diagram starting at the onset of 
mass transfer, for the Z = Zq Systems No. 1, 2, 3, 4, 7, 9, 
13 (identifiable by the corresponding initial main sequence 
star masses, where the main sequence star of two systems 
with the same initial main sequence star mass has a larger 
initial luminosity for the system with the larger initial 
period). The tracks are shown as continuous lines up to the 
point when the white dwarf reaches 1.4 Mq, and are then 
continued as dotted lines (cf. Sect. 3). The shaded band 
is limited by the dot-dashed lines connecting the zero age 
and the terminal age main sequence positions of single 
stars in the mass range O.8...3M0, taken from Schaller 
et al. (1992). The two dashed lines give the evolution of 
the nuclear luminosity Lnuc = Snucdm as function of 
the main sequence star mass from the onset to the mass 
transfer to the time when the white dwarf reaches 1.4 Mq, 
for Systems No. 1 (lower line) and No. 2 (upper line). 

these stars are very uncvolvcd, i.e. their internal hydrogen 
and helium distribution is basically homogeneous. There- 
fore, when these stars have lost a significant part of their 
initial mass, the spatial distribution of their main con- 
stituents (H and He) is still the same as in a zero age main 
sequence star. Since the stellar structure adjusts according 
to the chemical profiles, and once the mass transfer rate 
has dropped and the star can relax to thermal equilib- 
rium, their properties can not be distinguished from those 
of normal main sequence stars of the same mass (except 
for trace elements; see below). 

The deviation of these stars from the zero age main 
scqiience mass-himinosity relation is due to the thermal 
imbalance. Due to the strong mass loss, the outer stellar 
layers expand which consumes energy and reduces the stel- 
lar luminosity as long as the strong mass loss prevails — 
strong meaning that the mass loss time scale is of the same 
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Fig. 19. Same as Fig. 18, but for our Z = 0.001.^0 Sys- 
tems No. 60, 61, 63, 64, 65, 67, 68, and 70. The evolution of 
the nuclear luminosity is shown for Systems No. 60 (lower 
line) and No. 61 (upper line). 

order as the stellar Kelvin-Helmholtz time scale. This ef- 
fect by itself can reduce the stellar luminosity by as much 
as a factor of 10 (Figs. 18 and 19), i.e., the star can appear 
ten times dimmer than expected from the mass-luminosity 
relation of single stars on the basis of its actual mass. The 
effect is stronger for larger mass loss rates. We empha- 
sise that this reduction of the stellar luminosity due to 
the mass loss induced thermal imbalance comes actually 
in two components. One is that due to the strong mass 
loss, the outer stellar layers expand which consumes en- 
ergy and reduces the stellar luminosity below the "nuclear 
luminosity" , i.e. the amount of energy liberated by ther- 
monuclear reactions in the stellar core. However, as main 
sequence stars adjust their nuclear luminosity to the ra- 
diative energy loss at the surface (Kippenhahn & Weigert 
1990), also the nuclear luminosity of the mass losing main 
sequence stars is smaller than the nuclear luminosity of a 
non-mass losing main sequence star of the same mass and 
evolutionary stage (see Figs. 18 and 19). 

For systems with a relatively large initial period, i.e. 
those where the mass transfer starts when the main se- 
quence star is close to the terminal age main sequence 
(upper borderline of the main sequence band in Figs. 18 
and 19), one phenomenon counterbalances the two dim- 
ming effects (i.e. mass reduction and thermal imbalance). 
The cores of these stars are, at the onset of the mass trans- 
fer, very helium-rich. Therefore, after the mass transfer 
they have a helium-rich core which is significantly more 
massive than a helium core in a single stars of the same 
stellar mass. This makes the stars overluminous compared 
to single stars, as can be seen from the dotted parts of the 
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Fig. 20. Evolutionary tracks of the main sequence stars 
in the mass-radius diagram starting at the onset of mass 
transfer, for the Z = Zq Systems No. 1, 2, 3, 7, 9 (iden- 
tifiable by the corresponding initial main sequence star 
masses, where the main sequence star of two systems with 
the same initial main sequence star mass has a larger ini- 
tial radius for the system with the larger initial period). 
The tracks are shown as continuous lines up to the point 
when the white dwarf reaches 1.4 Mq, and are then con- 
tinued as dotted lines (cf. Sect. 3). The shaded band is 
limited by the dot-dashed lines connecting the zero age 
and the terminal age main sequence positions of single 
stars in the mass range O.8...3M0, taken from Schaller et 
al. (1992). 



stellar tracks of the long-period systems shown in Figs. 18 
and 19. Evidently, this effect — we call it the helium effect 
— is the larger the later the mass transfer starts during 
the core hydrogen burning evolution of the main sequence 
component, and the larger the total amount of mass lost. 

In summary, we have three effects changing the lu- 
minosity of our main sequence component. First, as its 
mass is reduced, its luminosity is reduced according to 
the mass-luminosity relation of single stars (the mass ef- 
fect). Second, the larger the mass loss or mass transfer 
rate the more is its luminosity reduced in addition, due 
to the thermal imbalance imposed by the mass loss (the 
mass loss effect). Third, the helium effect can lead to an 
increase of the luminosity for those stars which started 
out with (and therefore still have) relatively large peri- 
ods. All together, we see from Figs. 18 and 19 that during 
the thermally unstable mass transfer phase no star can be 
found above the single star mass-luminosity band. On the 
other hand, a significant fraction of them is found below 
this band, at luminosities between 1 and 10 Lq. 
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Fig. 20 shows, for the case of the metal-rich stars, that 
not only the luminosities but also the radii of the mass 
transferring main sequence stars are significantly smaller 
than the radii of stars in thermal equilibrium with the 
same mass and evolutionary stage. Note that there is also 
a strong dependence of the main sequence star radii on 
metallicity (cf. Figs. 1 and 2 in Sect. 2.3). Both effects may 
be quite relevant for the derivation of component masses 
in supersoft X-ray binaries (cf. Sect. 5.1). 

5. Observable consequences 
5.1. Supersoft X-ray sources 

White dwarfs which accrete hydrogen at such a rate that 
they can perform non-explosive hydrogen burning at their 
surface constitute the leading model for the persistent su- 
persoft X-ray sources (SSSs, Kahabka & van den Heuvel 
1997). Here, we want to compare our results to observa- 
tions of SSSs. It is important to keep in mind that we 
restricted the parameter space of our models according 
to the possibility to obtain a Type la supernova. I.e., al- 
though all our models may be considered as models for 
SSSs, it is not excluded that the average SSSs have in 
fact quite different properties than our models. We there- 
fore restrict ourselves to investigate three basic observable 
properties, i.e. the X-ray luminosity, the system period, 
and the luminosity of the donor star, and rather focus on 
what the largest and smallest of these values are rather 
than considering a typical average. For this purpose, we 
have compiled in Table 4 system properties at the time of 
the maximum X-ray luminosity — i.e., at the time of the 
maximum mass accumulation rate of the white dwarf, as 
we assume Lx = eMwd (cf. Sect. 2.1). 

The maximum X-ray luminosity which we can achieve 
in principle within our assumptions is that obtained 
by a Chandrasekhar-mass white dwarf accreting at 
its Eddington-rate (cf. Sect. 2.1), i.e. 2.07 lO'"'^ erg s^^ 
The largest value actually occurring in our models is 
1.85 10"^* erg s^^ (cf. Table 4). So far, none of the empiri- 
cal bolometric fluxes derived from SSSs exceeds this value, 
although some are quite close to it (Kahabka & van den 
Heuvel 1997). 

For comparing the periods of our models with those 
of SSSs, we focus on the short periods, since also post 
main sequence donor stars can produce SSSs which would 
occur in longer period systems (Li & van den Heuvel 1997, 
Hachisu et al. 1999, Wellstein et al. 1999b). At Z = 0.02, 
we find periods in the range 1.8...0.25d (43...6h), while at 
Z = 0.001 periods range from 19...5h (cf. Tables 2 and 3). 
Observed periods in close binary SSSs (Kahabka & van 
den Heuvel 1997) are generally in good agreement with 
these figures. 

Some authors in the literature express the necessity to 

explain the shortest periods found in SSSs with alternative 
scenarios. E.g., the SMC system 1E0035.4-7230 has a pe- 



riod of 4.1 h (Schmidtke et al. 1996), for which van Teesel- 
ing & King (1998) proposed a wind-driven evolution, with 
a very low mass main sequence star losing mass induced 
by the strong X-ray radiation of the white dwarf. We note 
that in particular our low-Z models show periods as low as 
4h (e.g.. System No. 71 in Table 3). Furthermore, accord- 
ing to Eq. (19) significantly smaller periods are achievable 
for smaller but still plausible initial mass ratios. I.e., a sys- 
tem starting out with a 2.4 Mq main sequence star and a 
0.6 Mq white dwarf could reduce is initial period by a fac- 
tor 4. Thereby, even periods in the range 2...3h could be 
obtained. Even though such systems might not lead to 
Type la supernovae, some of them may still allow for sta- 
tionary hydrogen burning on the white dwarf surface for 
a limited amount of time. 

Rappaport et al. (1994), in a population synthesis 
study of SSSs, considered all possible initial masses and 
periods. The shortest periods they find are of the order of 
5h. They consider only one metallicity (solar), for which 
the smallest period we found is ~ 6h. In our low metal- 
licity systems, we find a minimum of 4 h. Thus, also from 
to the results of Rappaport et al. we would expect then 
minimum periods as low as 3 h at low metallicity. There- 
fore, periods as short as that found in 1E0035.4-7230 may 
still be explained within the standard model of thermally 
unstable mass transfer studied in the present paper. 

Finally, we want to discuss the brightness of the donor 
stars in our models, in relation to the fact that so far none 
of them could be observationally identified. In Sect. 4.2 
we have seen that the main sequence stars in our models 
are, during the mass transfer phase, significantly imder- 
luminous for their actual mass. In Table 4, we show the 
properties of the main sequence stars at the time of the 
maximum X-ray luminosity, for selected cases. Comparing 
Systems No. 1 and No. 3, we see that the stellar luminosity 
is not well correlated with the mass of the main sequence 
star. Instead, it is inversely correlated with the mass trans- 
fer rate, i.e. the mass loss rate of the main sequence star. 
For systems which have no wind, i.e. for which M = Mwd 
and thus Lx oc M, this means that the brighter the sys- 
tem in X-rays, the dimmer is the main sequence star. I.e., 
the fact that in the supersoft X-ray sources the X-ray lu- 
minosity is large (otherwise we would not notice them) 
means that the mass transfer rate must also be large (cf. 
Sect. 4.3). We conclude that the reduction of the main 
sequence star luminosity due to the thermal imbalance 
must be a large effect in observed supersoft sources of the 
considered type. As it can reduce the bolometric luminos- 
ity of the main sequence star by more than one order of 
magnitude, it may be quite difficult to observe the main 
sequence component in supersoft X-ray binaries. 

5.2. The stellar remnant 

Once the white dwarf has exploded, the main sequence 
component is likely to survive, and although some small 
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Table 4. Properties of selected systems at the time of the maxinmm X-ray luminosity. The columns have the following 
meanings. (1) system number (cf. Tables 2 and 3), (2) main sequence star initial mass, (3) white dwarf initial mass, 
(4) initial period, (5) mass transfer rate, (6) mass accumulation rate of white dwarf, (7) white dwarf luminosity (8) 
Eddington luminosity of white dwarf (9) effective temperature of white dwarf, (10) white dwarf mass, (11) main 
sequence star mass, (12) main sequence star luminosity, (13) effective temperature of main sequence star, (14) orbital 
period (15) orbital velocity of main sequence star, (16) orbital velocity of white dwarf. 



Nr. 


Mms.i 


MwD,i 


Pi 


M 


MwD 


Lx 


^Edd 


2wD 


MwD 


Mms 


Lms 


Tms 


P 


VMS 


VWD 










lo-'^ 


lo-'^ 


1q38 


1q38 


10^ 








10^ 










Mq 


M© 


d 


Mo/yr 


Mo/yr 


erg/s 


erg/s 


K 


Me 


Mq 


L© 


K 


d 


km/s 


km/s 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


1 


2.3 


1.0 


0.51 


7.47 


3.87 


1.57 


1.57 


818 


1.06 


1.91 


2.00 


6.26 


0.32 


160 


288 


3 


2.1 


1.0 


1.65 


3.90 


3.90 


1.85 


1.85 


877 


1.25 


1.83 


10.7 


6.17 


1.23 


117 


172 


4 


2.0 


1.0 


0.69 


3.75 


3.75 


1.66 


1.66 


837 


1.12 


1.87 


7.03 


7.02 


0.59 


137 


229 


9 


1.8 


1.0 


0.54 


1.45 


1.45 


0.64 


1.64 


648 


1.11 


1.67 


5.01 


7.09 


0.47 


154 


231 


29 


2.1 


0.8 


1.19 


2.39 


2.39 


1.05 


1.51 


734 


1.02 


1.43 


1.90 


5.73 


1.14 


114 


160 


34 


2.0 


0.8 


1.55 


2.38 


2.38 


1.05 


1.57 


739 


1.06 


1.60 


5.56 


5.96 


1.37 


106 


160 


36 


1.9 


0.8 


0.60 


2.40 


2.38 


1.05 


1.30 


744 


0.88 


1.82 


6.35 


7.10 


0.54 


118 


245 


38 


1.8 


0.8 


0.58 


2.08 


2.08 


0.92 


1.30 


756 


0.88 


1.57 


5.37 


7.01 


0.50 


130 


232 


40 


1.7 


0.8 


0.66 


1.40 


1.40 


0.62 


1.41 


636 


0.95 


1.54 


3.98 


6.58 


0.50 


139 


225 


49 


2.0 


0.7 


0.47 


8.05 


2.67 


1.20 


1.20 


731 


0.81 


1.71 


1.26 


6.07 


0.25 


148 


312 


51 


1.9 


0.7 


0.48 


8.60 


2.87 


1.28 


1.28 


751 


0.87 


1.46 


0.86 


5.73 


0.25 


167 


280 


54 


1.9 


1.0 


0.31 


6.92 


3.70 


1.64 


1.64 


833 


1.11 


1.64 


1.80 


6.90 


0.23 


196 


290 


60 


1.8 


1.0 


0.29 


5.21 


3.81 


1.69 


1.69 


843 


1.14 


1.60 


2.17 


7.27 


0.24 


200 


280 


62 


1.7 


1.0 


0.29 


4.07 


3.86 


1.70 


1.70 


845 


1.15 


1.54 


2.58 


7.54 


0.25 


127 


170 


64 


1.6 


1.0 


0.29 


2.85 


2.85 


1.25 


1.69 


783 


1.14 


1.46 


2.65 


7.62 


0.25 


204 


261 


67 


1.5 


1.0 


0.61 


3.36 


3.36 


1.48 


1.70 


817 


1.15 


1.35 


5.62 


7.25 


0.54 


164 


192 


70 


14 


1.0 


0.58 


2.09 


2.09 


0.92 


1.67 


724 


1.13 


1.26 


5.04 


7.19 


0.53 


167 


186 


71 


1.8 


0.8 


0.24 


9.73 


2.38 


1.06 


1.36 


725 


0.92 


1.50 


0.65 


5.97 


0.18 


194 


314 



amounts of mass may be stripped off by the supernova 
ejecta and blast wave, (Wheeler et al. 1975, sec also Fryx- 
ell & Arnett 1981, Taam and Pryxell 1984), most stellar 
properties of our main sequence components will remain 
more or less unchanged. One may hope to identify and ob- 
serve the remaining main sequence component either in a 
young Galactic supernova remnant produced by a Type la 
supernova, or, if they stick out sufHciently, long after the 
supernova explosion in the field. 

In the first case, the thermal imbalance imposed by 
the mass transfer (cf. Sect. 4.4) will still be completely 
preserved, since the thermal time scale of the star is of 
the order of 10^ yr, while any gaseous supernova remnant 
would dissolve at least 100 times faster. The expected lu- 
minosities of the main sequence components in a super- 
nova la remnant can thus be directly read off Figs. 18 
and 19 from the star-symbols, which mark the expected 
time of the supernova explosion. They are found to be in 
the range 1...10Lq. When inspecting the effective tem- 
peratures of the main sequence stars at the time of the 
supernova, we find them to be systematically 500...1000K 
cooler than single main sequence stars of the same mass 
and evolutionary state (Fig. 17; see also column 14 in Ta- 
ble 4). As the main sequence band has a width of more 
than 1000 K, this implies that remnant stars will be lo- 
cated on the main sequence band or slightly to the right. 



At Z = 0.02, the effective temperatures of the remnant 
stars are larger than 5500 K, at Z = 0.001 larger than 
6000 K. I.e., they would appear as evolved F or G type 
main sequence stars. 

Important to unambiguously identify the stellar rem- 
nant of a supernova la progenitor system of the consid- 
ered type is its peculiar surface chemical composition. 
These stars have peculiar abundances since they have 
lost a major part of their initial mass during the mass 
transfer phase, with the consequence that they uncover 
matter which has been sufficiently deep inside the star 
that thermonuclear reactions have occurred.. All main se- 
quence stars in the present study in systems which lead to 
Type la supernovae lose at least ^ 0.4 M0, as the white 
dwarf needs to achieve the Chandrasekhar mass. However, 
in those systems where the white dwarf develops a wind 
the total mass loss of the main sequence stars may be 
considerably larger (cf. Table 5). 

In a normal main sequence star, all isotopes of the 
light elements lithium, beryllium and boron are destroyed 
in the whole stellar interior except in an outer envelope of 
~ 0.1 Mq (or ~ 0.2 Mq for boron). Therefore, the main 
sequence components of our systems are, at the time the 
supernova explosion occurs, all completely devoid of the 
light elements. The lack of these elements offers there- 
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Table 5. Ratios of surface abundances of the main sequence star to initial abundance, for the time of the supernova 
explosion. The abundances of all isotopes of the light elements LI, Be, and B are zero for all models. 



Nr. 


MMS.i 

Mo 


MMS,i 
Mq 


AM 

Mq 


AMwind 


^He 


^He 


12^ 


13^ 




15N 


16q 


"0 


18q 


^^Na 


3 


2.1 


1.66 


0.44 


0.00 


6.71 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


8 


2.0 


1.55 


0.45 


0.00 


7.91 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


13 


1.8 


1.35 


0.45 


0.00 


14.8 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


29 


2.1 


0.74 


1.36 


0.72 


4.33 


1.00 


0.08 


2.22 


4.61 


0.06 


1.00 


2.92 


0.02 


1.06 


31 


2.1 


0.79 


1.31 


0.66 


6.65 


1.04 


0.15 


3.66 


4.28 


0.07 


1.00 


2.39 


0.05 


1.07 


32 


2.0 


1.05 


0.95 


0.31 


7.45 


1.01 


0.91 


7.28 


1.12 


0.34 


1.00 


1.00 


0.77 


1.00 


35 


2.0 


0.66 


1.34 


0.69 


4.10 


1.07 


0.08 


1.91 


4.60 


0.06 


0.99 


6.27 


0.03 


1.07 


36 


1.9 


1.17 


0.73 


0.09 


16.8 


1.01 


0.99 


1.48 


1.00 


0.89 


1.00 


1.00 


0.99 


1.00 


37 


1.9 


0.93 


0.97 


0.33 


14.4 


1.03 


0.69 


10.7 


1.80 


0.22 


1.00 


1.02 


0.46 


1.00 


38 


1.8 


1.15 


0.65 


0.00 


22.6 


1.00 


0.99 


1.10 


1.00 


0.97 


1.00 


1.00 


0.99 


1.00 


39 


1.8 


1.15 


0.65 


0.00 


26.9 


1.01 


0.99 


1.30 


1.00 


0.93 


1.00 


1.00 


0.99 


1.00 


45 


1.7 


1.06 


0.64 


0.00 


28.1 


1.01 


0.99 


1.64 


1.00 


0.90 


1.00 


1.00 


0.99 


1.00 


54 


1.9 


1.25 


0.65 


0.21 


218 


1.00 


0.98 


2.42 


1.00 


0.71 


1.00 


1.00 


0.96 


1.00 


55 


1.9 


1.12 


0.78 


0.34 


172 


1.02 


0.80 


11.1 


1.38 


0.20 


1.00 


1.01 


0.55 


1.00 


60 


1.8 


1.31 


0.49 


0.04 


209 


1.00 


0.99 


1.05 


1.00 


0.98 


1.00 


1.00 


1.00 


1.00 


61 


1.8 


1.07 


0.73 


0.29 


281 


1.02 


0.85 


9.80 


1.26 


0.24 


1.00 


1.01 


0.65 


1.00 


62 


1.7 


1.25 


0.45 


0.00 


195 


1.00 


0.99 


1.02 


1.00 


0.99 


1.00 


1.00 


1.00 


1.00 


63 


1.7 


1.23 


0.47 


0.13 


474 


1.00 


0.97 


1.98 


1.00 


0.79 


1.00 


1.00 


0.97 


1.00 


64 


1.6 


1.15 


0.45 


0.00 


208 


1.00 


0.99 


1.02 


1.00 


0.99 


1.00 


1.00 


1.00 


1.00 


65 


1.6 


1.13 


0.47 


0.03 


622 


1.00 


0.99 


1.13 


1.00 


0.97 


1.00 


1.00 


0.99 


1.00 


67 


1.5 


1.06 


0.44 


0.00 


706 


1.00 


0.99 


1.09 


1.00 


0.98 


1.00 


1.00 


1.00 


1.00 



fore already an unambiguous way to identify the remnant 
stars. 

In Table 5, we compile other surface abundance 
anomalies found in our models. It can be seen, that the 
isotope ^He is overabundant by a large factor, which is, 
however, hard if not impossible to verify observationally 
at the present time. The same may hold for other isotopic 
anomalies, e.g. of ^^C and ^^N. Only for carbon and ni- 
trogen, we find the possibility of peculiar elemental abun- 
dances, i.e., carbon may be significantly underabundant 
and nitrogen correspondingly overabundant, as the CN- 
cycle is responsible for this feature. 

The supernova explosion can in principle alter the 
abundances shown in Table 5 in two ways. It can lead 
to an additional mass loss of the main sequence star of 
the order of 0.1 Mq...0.2 Mg (Wheeler et al. 1975, Fryxell 
& Arnett 1981, Taam and Fryxell 1984). This may lead 
to a somewhat stronger CN-cycle signature but would not 
change our results qualitatively. It is further not excluded 
that it can lead to the deposition of small amounts of the 
supernova ejecta on the main sequence star (Fryxell & 
Arnett 1981, Taam and Fryxell 1984). Whether this hap- 
pens or not seems to be unclear at the moment. In any 
case, the effect might be some enrichment of the surface 
composition of the main sequence star remnant with the 
nucleosynthesis products of the supernova, i.e. in elements 
between carbon and iron (Thielemann et al. 1986). 

In summary, the light elements, e.g., lithium, and 
carbon are the most promising distinguishing chemical 



characteristics of main sequence type stellar remnants of 
Type la supernovae. Another independent characteristic 
may be a peculiar radial velocity or proper motion. The 
main sequence stellar remnants will at least have a pe- 
culiar velocity of the order of their orbital velocity at the 
time of the explosion of the white dwarf. This velocity is in 
the range 140kms~^...250kms~^ for all our systems, with 
larger values corresponding to the low metallicity models. 
The momentum impacted by the supernova ejecta on the 
main sequence star may increase its space velocity up to 
as much as ~ 500kms~^ (Wheeler et al. 1975). Therefore, 
any main sequence type remnant star must have space ve- 
locities in a very favourable range. It is large enough to 
impose a clearly peculiar kinematic on the stellar rem- 
nant, but it is still much smaller than the velocity of the 
supernova ejecta, which implies that the star will remain 
for a long time close to the center of the supernova rem- 
nant. Chemical and kinematic signature together make it 
in fact a interesting project to search for a main sequence 
type stellar remnant in the gaseous remnant of the his- 
torical Galactic supernova 1006 (Wellstein et al. 1999b), 
which is very likely the product of an exploding white 
dwarf (Schaefer 1996). 

6. Discussion and conclusions 

We have studied the evolution of close binary systems con- 
sisting of a main sequence star and a white dwarf which 
are considered as candidates for progenitors of Type la su- 
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pcrnovac. Based on an extended grid of models, we have 
studied the properties of the systems as a function of the 
initial donor star mass, initial white dwarf mass, initial 
period, and chemical composition. 

Due to our numerical technique (Sect. 2) we obtain, for 
the first time, a complete picture of the time dependence 
of the mass transfer rate in such systems. We find that the 
mass transfer rate remains initially for about one thermal 
time scale of the main sequence star on a very low level 
during which nova outbursts are likely to occur. Then, the 
maximum mass transfer rate is rapidly reached. We find 
that most white dwarfs approach the Chandrasekhar mass 
during the decline phase of the mass transfer (A^wd < 0; 
cf. Figs. 7 to 10). Our results will allow to investigate the 
effect of this time dependence of the mass transfer rate 
on the lower critical accretion rates for stationary nuclear 
burning on the white dwarf (cf. Sect. 2.1). This may be 
important, as Prialnik & Kovetz (1995) showed that these 
threshold values may be smaller for higher white dwarf 
temperatures. As for Mwd < the white dwarf temper- 
ature is expected to be higher at a given value of Mwd 
compared to the case of Mwd = 0, this effect may per- 
haps increase the parameter space of models which lead 
to Chandrasekhar mass white dwarfs. 

In contrast to results based on simple estimates of 
the mass transfer rate (e.g., Eq. (13)), we find that the 
rates increase strongly for lower initial white dwarf masses 
(Fig. 9). I.e., even systems with rather small initial white 
dwarf masses (~ 0.7 M©) can not be excluded to evolve to 
Type la supernovae. As Chandrasekhar mass white dwarfs 
are likely to rotate faster the smaller their initial mass is, 
this implies that the white dwarf rotation may be relevant 
in Type la explosions (cf. Sect. 4.1). 

We find that the mass transfer rates in low metallicity 
systems are, for the same initial main sequence star and 
white dwarf masses, much higher than at solar metallic- 
ity. I.e., the initial main sequence star mass range which 
results in good Type la supernova progenitor candidates 
shifts from 1.6 M©... 2.3 M© to I.4M0...I.8M0 (Figs. 14 
and 15). We note that the exact donor star mass range is 
uncertain, due to uncertainties in white dwarf wind mass 
accumulation efficiencies (cf. Sects. 2 and 3.2). However, 
we find that at low metallicity, this range is narrower, 
and supernova la progenitor systems need to have white 
dwarfs which are initially about 0.2 Mq more massive (see 
also Figs. 14 and 15), leading to a decrease of the Type la 
supernova rate with decreasing metallicity. We note that 
this effect differs from the wind inhibition effect proposed 
by Kobayashi et al. (1998). 

It is of course tempting to speculate about effects of 
the metallicity dependence of the progenitor evolution on 
the supernova peak brightness or decline rate. However, 
although more realistic calculations of the white dwarf 
evolution to the Chandrasekhar mass are now possible, 
they need to be performed before definite conclusions can 
be drawn. This is so since, although the described effects 



are likely to introduce a Z-depcndence to the supernova 
properties, it would interfere with other such effects as 
described by Hoflich et al. (1998), Dominguez et al. (1999), 
and Umcda ct al. (1999). 

It is also hard to disentangle whether possible depen- 
dences of Type la supernova properties on their environ- 
ment — as suggested by Branch et al. (1996) or Wang et 
al. (1997) — are due to the mentioned trends with metal- 
licity or due to different progenitor types at work. The 
latter seems more likely considering the life times of our 
progenitor models (Tables 2 and 3). Although our low Z 
models invoke lower mass main sequence stars, their pre- 
dicted life time is < 1.5 10® yr, which is roughly similar for 
our solar metallicity models. I.e., Type la supernovae in 
elliptical galaxies, which may require progenitor life times 
of 10"'^° yr, can not be obtained from the type of model 
presented here, but may rather require scenarios with low 
mass red giant (Hachisu et al. 1996) or CO white dwarf 
(Iben & Tutukov 1984) donor stars. 

We emphasise that, nevertheless, models of the consid- 
ered type very likely exist in nature, as they correspond 
to the close binary supersoft X-ray soTirces (Kahabka & 
van den Heuvel 1997). The X-ray luminosities, periods, 
and main sequence star properties (Sect. 5.1 and Table 4) 
appear to agree quite well with observed systems. We also 
outline a way to test whether such systems can in fact 
evolve into Type la supernovae. We make unambiguous 
predictions for the chemical and kinematical properties of 
the stellar remnants of main sequence star + white dwarf 
systems after the explosion of the white dwarf (Sect. 5.2), 
which may be directly tested for the case of the historical 
galactic supernova 1006 (Wellstein ct al. 1999b). 
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